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Escape of Lya and continuum photons from star-forming 
galaxies 
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ABSTRACT 

A large number of high-rcdshift galaxies have been discovered through narrow- 
band Lya line or broad-band continuum in recent years. The escaping process of 
photons from these early galaxies is crucial to understanding galaxy evolution and 
the cosmic rcionization. Here, we investigate the escape of Lya, non-ionizing UV- 
continuum (A = 1300 — 1600 A in rest frame), and ionizing photons (A < 912 A) from 
galaxies by combining cosmological hydrodynamic simulations with three-dimensional 
multi-wavelength radiative transfer calculations. We find that the escape fraction (/ esc ) 
of these different photons shows a complex dependance on redshift and galaxy prop- 
erties: f^. a and f^f appear to evolve with redshift, and they show similar, weak 
correlations with galaxy properties such as mass, star formation, mctallicity, and dust 
content, while /*°" remains roughly constant at ~ 0.2 from z ~ — 10, and it does 
not show clear dependence on the properties of the galaxy. The f^/ c a correlates more 
strongly with /^J than with J'™. In addition, we estimate the ionizing photon emis- 
sivity of Lyman Alpha Emitters (LAEs) and their contribution to the ionization of 
intcrgalactic medium (IGM), by combining our simulations with the observed luminos- 
ity functions of LAEs at different redshift. Our result suggests that ionizing photons 
from LAEs alone are not sufficient to ionize IGM at z > 6, but they can maintain the 
ionization of IGM at z ~ — 5. 



Key words: radiative transfer - ISM: dust, extinction - galaxies: evolution - galaxies: 
formation - galaxies: high-redshift 



1 INTRODUCTION 

Young, star forming galaxies can produce strong Lya 
emission via hydrogen reco mbination in the ionized region 
|Partridge fc Peeblesl 1 19671 ). A large population of Lya 
emitters, or LAEs, has been detected over a wide redshift 
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Therefore, Lya emission may be 
a powerful tool to study galaxy evolution. 

The intrinsic Lya flux is proportional to star formation 
rate (SFR) when the ionizing photons are absorbed by neu- 
tral hydrogen. However, most Lya photons can experience 
numerous scattering processes by neutral hydrogen due to 
its large cross section. As a result, Lya photons can have a 
large optical depth and can be efficiently absorbed by inter- 
stellar dust. Hence, although most star forming galaxies are 
LAEs intrinsically, some may not show a strong Lya line due 
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to efficient dust absorption. It is therefore crucial to study 
the escape of Lya photons to understand LAEs. 

Observationally, some groups have derived the escape 
fraction of Lya from th e observed flux ratio be- 

tween Ly a and Ha at z = 0.3 l|Atek et al.l 120091 ) . and at 
z = 2.2 (|Haves et alj|201Ct ). It was shown that / e L s ? Q de- 
creases with color excess _E(B — V), albeit with a large dis- 
persion (|Atek et al.l 120091 ) . Alternatively , f^ r a can be de- 
rived from Lya a nd UV continuum flux (|Haves et alj|201ll ; 
iBlanc et ailkoill ). and it was suggested that it monotoni- 
cally increases with redshift and can be fitted with simple 
power-law function (|Haves et alj|201lT h 

Theoreticall y, semi-analytical studies (e.g., 

iKobavashi et al] l2010h and cosmological h ydrodynam- 
Nagamine et al.l |2010| ; iDaval et al.l 



ics simulations (e 



120091 ; [Sh muzu et alll201ll ) indicate that an adjustment of 
Lya flux by a factor of / csc is needed to fit the observed 
luminosity function or Lya equivalent width distribution. 
More recently, by combining cosmological hydrodynamics 
and three-dimensional Lya radiative transfer (RT) simu- 
lations, it was suggested that the es cape fraction of Lya 
photons decreases with halo mass (|Laursen et al.l 120091 ; 
lYaiima et al1l201lh . 

Moreover, the ionizing photons escaping from LAEs 
may contribute significantly to the cosmic reionization. It 
is widely believed that the inter-galactic medium (IGM) 
is highly ionized z < 6 from the Gunn- Peterson troughs 
(|Gunn fc Peterson] 1965) in the observed quasar (QSO) 
spectra fe.g.. iFan et al] 2008) , and that it should be par- 



tially ionized even at z > 10 from the Thomson scattering 
optic al depth of the cosm ic microwave background radia- 
tion (|Komatsu et al]|2009l ). The possible ionizing sources 
are thought to be star-forming galaxies, QSOs, and PopIII 
stars. However, since the observed number density of QSOs 
drops drastically at z > 4, the ionizing photons from QSOs 
alone ca nnot maintain the ioniza tion of the IGM at high 
redshift (|Bolton fc Haehnerj2007l ) . On the other hand, since 
the number density of LAEs d oes not decrease significantly 
at 2 > 3 (|Ouchi et al] |2008l ) , and the contribution from 
LAEs to the cosmic star formation rate densit y becomes 
dominant at high redshift (|Ciardullo et al]|2012h . so high-z 
LAEs may play an important role in the ionization of IGM. 

However, the estimation of ionizing ability of LAEs in 
previous work suffered from the ambiguities in the escape 
fraction of Lya and ionizing photons. Moreover, an equiva- 
lent width (EW) of larger than 20 A is frequent ly used as 
the definition of LAEs (e.g- jGronwall et alj|2007h . Hence, it 
is important to understand the difference of escape fraction 
between Lya and non-ionizing UV continuum. As the pho- 
ton escaping depends sensitively on the galaxy properties, 
the ionization state of the interstellar medium (ISM), and 
the dust content, it is necessary to investigate the radiative 
transfer of these different photons in galaxies simultaneously. 

In this work, we study the escaping of Lya, non-ionizing 
UV-continuum (A = 1300 — 1600 A in rest frame, hereafter 
simply "UV continuum"), and ionizing photons (A < 912 A) 
from galaxies by combining cosmological hydrodynamic sim- 
ulations with three-dimensional multi-wavelength radiative 
transfer calculations. The paper is organized as follows. We 
describe the galaxy model in §2, and our multi-wavelength 
radiative transfer code ART 2 in §3. In §4, we present the re- 
sults of the escape of Lya, UV-continuum, and ionizing pho- 



tons, including their evolution with redshift, dependence on 
galaxy properties, and the contribution from LAEs to cos- 
mic reionization. We discuss the dependence of / osc on the 
dumpiness of the ISM, and the limitations of our calcula- 
tions in §5, and summarize our finding in §6. 



2 GALAXY MODEL 

The cosmological simulation presented here follows the for- 
mation and evolution of a Milky Way-like galaxy and its sub- 
structures. The si mulation has been emp loyed in a series of 
RT calc ulations bv lYaiima et al] (|2012al ) and lYaiima et ail 
(|2012bl ). and it will be described in detail in Zhu et al. (in 
preparation). Here we briefly highlight some important fea- 
tures of the simulation. 

It uses the Aquila initial condition, which is the hy- 
drodynamical version (jWadepuhl fe Springell l201ll ) of the 
initial condition of the Aquarius Project, the largest ever 
particle simulation o f a Milky Way-sized dark matter halo 
(|Springel et al] 120081 ) . The simulation includes dark mat- 
ter, gasdynamics, star formation, black hole growth, and 
feedback processes. It was performed using the paral- 
lel, N-body/Smoothed Particle Hydrodynamics (SPH) code 
GADGE T-3, which is a n im prove d version of t hat de- 
scribed in lSpringel et al] (|200ll ) and lSpringeil (|2005l ). GAD- 
GET implements the entropy-co nserving formulation of 
SPH (jSpringel fc Hernquistl I2002T) with adaptive particle 
smoothing, as in iHernquist fc Katzl (|l989h . Radiative cool- 
ing and heating processes a re calculated assuming colli- 
siona l ionization equilibrium (|Katz et al] 1 19961 : iDave et al] 
Il999h . Star formation is modeled in a multi-phase I SM, with 
a rate that follows th e Schmidt-Kennicutt Law (jSchmidtl 
Il959l : lKennicutdll998h . Feedback from supernovae is cap- 
tured through a multi-phase model of t he ISM by an effective 
equa tion of state for star-forming gas (ISpringel fc Hernquist 
120031) . The UV background model of lHaardt fc Madau 
(|l996h is used. 

The whole simulation falls in a periodic box of 
100 /i~ 1 Mpc on each side with a zoom-in region of a size 
5x5x5 /i~ 3 Mpc 3 . The spatial resolution in the zoom-in re- 
gion is ~ 250 h^ 1 pc in comoving scale. The mass resolution 



of this zoom-in region is 1.8 x 10 B h 1 Mq for dark mat- 
ter particles, 3 x 10 5 h^M® for gas, and 1.5 x 10 5 h _1 M Q 
for star particles. The cosmological parameters used in the 
simulation are Q m = 0.25, Qa — 0.75, ag = 0.9 and 
h = 0.73, consistent w ith the five-year results of the WMAP 
i|Komatsu et al]|2009l ') . The simulation evolves from z = 127 
to z = 0. 

In the simulation, a galaxy is identified as a group using 
the on-fly friends-of-friends (FOF) group finding algorithms, 
which links all baryon particles to their nearest dark matter 
neighbor with a dark matter linking length less than 20% of 
their mean spacing. 



3 THE RADIATIVE TRANSFER CODE ART 

The RT calculations are performed using the 3D Monte 
Carlo RT code, All-wavelength Radiative Trans fer with 
Adaptive Refinement Tree (ART 2 ), as described in lLi et all 
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l|200Sh and lYaiima et all (|2012al ). ART 2 features three es- 
sential modules: continuum emission from X-ray to radio, 
Lya emission from both recombination and collisional exci- 
tation, and ionization of neutral hydrogen. The coupling of 
these three modules, together with an adaptive refinement 
grid, enables a self-consistent and accurate calculation of the 
Lya properties, which depend strongly on the UV contin- 
uum, ionization structure, and dust content of the object. 
Moreover, it efficiently produces multi-wavelength proper- 
ties, such as the spectral energy distribution and images, 
for direct comparison with multi-band observations. The de- 
tailed impleme ntati ons of the ART 2 cod e are described in 
iLi et al.l (|2008l ) and lYaiima et al.l |2012al ). Here we focus on 
the Lya calculations and briefly outline the process. 

The Lya emission is generated by two major mecha- 
nisms: recombination of ionizing photons and collisional ex- 
citation of hydrogen gas. In the recombination process, we 
consider ionization of neutral hydrogen by ionizing radiation 
from stars, active galactic nucleus (AGN), and UV back- 
ground (UVB), as well as by collisions by high-temperature 
gas. The ionized hydrogen atoms then recombine and create 
Lya photons via the state transition 2P — > IS. The Lya 
emissivity from the recombination is 



f a aBhu a n e nmi, 



(1) 



where ob is the case B recombination coefficient, and 
f a is the average number of Lya photons produced 
per case B recombi nation. Here we use ob derived in 
iHui fc Gnedlr] (| 19971 ). Since the temperature dependence 
of fg is not strong, f a = 0.68 is assumed everywhere 



(IQsterbrock fc Ferlan o feoOrJ ). The product hv a is the en- 
ergy of a Lya photon, 10.2 eV. 

In the process of collisional excitation, high temperature 
electrons can excite the quantum state of hydrogen gas by 
the collision. Due to the large Einstein A coefficient, the 
hydrogen gas can occur de-excitation with the Lya emission. 
The Lya emissivity by the collisional excitation is estimated 
by 



COII 



(2) 



where Chyc 



Cl y 



the collisional excitation coefficient, 
10~ 17 exp(-fa/ a /fcT)r~ 1/2 ergs s" 1 cm 3 



jOsterbrock fc Ferlandll2006l ) 

Once the ionization structure has been determined, we 
estimate the intrinsic Lya emissivity in each cell by the sum 
of above Lya emissivity, e a = £™ + e^ ' 1 - 

In RT calculations, dust extinction from the ISM is 
included. The dust content is estimated according to the 
gas content and metallicity in each cell, which are taken 
from the hydrodynamic simulation. The dust-to-gas ratio 
of the MW is used where the metallicity is of Solar abun- 
dance, and it is linearly interpolate d for other metallic i ty. We 
adopt the dust size distribution of iTodini fc Ferraral (|200ll ) 
for solar metallicity and a M = 22 Mq SN model, as in 
Figure 5 in their paper. The size distribution is then com- 
bi ned with the dust absorption and scattering cross section 
of IWeingartner fc Drain! (|200ll ) to calculate dust absorp- 
tion opacity curves. We use the stellar populati on synthesis 
model of GALAXEV (|Bruzual fc Charlot|[200l ') to produce 
intrinsic spectral energy distributions (SEDs) of stars for 
a grid of metallicity and age, and we use a simple, broken 



power law for the AGN (|Li et al.l 12003 ) . A lSalpeterl (|l955l ) 
initial mass function is used in our calculations. 

In this work, we apply ART 2 to the selected galaxies 
from the cosmological simulation. In our post-processing 
procedure, we first calculate the RT of ionizing photons 
(A ^ 912 A) and estimate the ionization fraction of the 
ISM. The resulting ionization structure is then used to run 
the Lya RT to derive the emissivity, followed by the calcula- 
tion of non-ionizing continuum photons (A > 912 A) in each 
cell. Our fiducial run is done with N p \-, = 10 5 photon packets 
for each ionizing, Lya, and non-ionizing comp onents, which 
was de monstrated to show good convergence (|Yaiima et al.1 
l2012al fbT). The highest refinement of the adaptive grid corre- 
sponds to a cell size comparable to the spatial resolution of 
250 pc in comoving coordinate of the hydrodynamic simula- 
tion. The physical values in each cell are interpolated from 
neighbor gas particles with the weight of the Kernel func- 
tion. The mean uniform density in each cell is used for ART 2 
calculations in our models. 



4 RESULTS 

4.1 Evolution of the escape fraction of Lya and 
continuum photons 

The time evolution of the escape fraction of Lya, UV con- 
tinuum at 1300 ^ A ^ 1600 A in rest frame, and ionizing 
photons is shown in Figure [T] We find that the / csc of Lya 
(/esc") moderately evolves with redshift, and the median val- 
ues range from ~ 0.2 to 0.9. The feZ a somewhat de- 
creases with increasing redshift at z < 3, and then increases 
with redshift z > 3. At lower redshfit, since the amount of 
the gas and dust decreases due to star formation, so some 
fraction of the Lya photons can escape. 

The median £° remains nearly constant, ~ 0.3 at 
2 ;$ z ;$ 4, which is similar to the observed value and 
trend with redshift in L AEs from the HETDEX pilot 
survey (|Blanc et al.l l201ll ) . The large di spersion in / P L ,'° , 
which is similar to the observations (e.g.. lAtek et al.ll2009l : 
lHaves et al1l2010l : iBlanc et al]|201ll ) , may be caused by dif- 
ferent galaxy properties such as mass, metallicity and star 
formation, as we will discuss later in the paper. At z > 3, 
our result i s in br oad agreement with the implication of 
lHaves et all (|201ll ). although the escape fractions of our 
sample increase more slowly with redshift. 

At z < 3, however , our r esults differ from the empirical 
model of lHaves et~aH |201ll ), in which / e L s ^ is parameter- 
ized as a power-law function of redshift, /<^c Q K (1 + z ) k i 
with the best fitting power-law index k = 2.6. This sim- 
ple relation implicitly comes from the redshift evolution of 
metallicity i n galaxies suggested b y some theoretical pre- 
dictions (e.g, Kobavashi et al.|[2007h and observations (e.g., 
iProchaska et al.ll2003l) .~ The dust content can increase with 
metal in gas, and can effectively absorb Lya photons. On the 
other hand, the distribution of gas and dust can strongly af- 
fect the f^ c a (e.g., Yajima et al. 2012c). At low redshift, 
the mean gas density decreases, and a large fraction of gas 
is consumed by star formation, so f^Z" of our model galaxies 
does not decreases with redshift at z < 3, despite the higher 
metallicity compared to that in higher redshift. In observa 
tions, there is a large dispersion in the estimation of /, 



Lya 
esc 



© 2008 RAS, MNRAS OOO.rflfTTl 



4 Yajima et al 




Figure 1. Escape fraction of Lya (upper panel), UV contin- 
uum at 1300 ^ A ^ 1600 A in rest frame (middle panel), and 
ionizing photons (lower panel) as a function of rcdshift, in com- 
parison with estimates from observations. The black open circles 
represent individual galaxies from the simulated sample, red filled 
circles are the corresponding median values at each redshift, with 
the error bars showing the quartiles. The blue symbols and ar- 
rows represent the escape fraction derived from observations. The 
triangle, cross an d square symbols i n th e upper panel are from 
lAtek et al] j2009h . lHaves et al]j201Ch and lBlanc et a.1.1 l|201lf ). re- 
spectively. The triangl e, square, cross and c i rcle s y mbols in the 



middle panel are from 



middl e p. 
1120041 1, [l 



Bouwens et al 



Adelberger fc Steidell (|200( 
feOOQl) and lHathi et all I 



1. lOuchi et all 



2008 



respec- 
tively, which are derived by co mbining the observe d E(B — V) 
with Calzctti's extinction law jCalzetti et al.ll2000h . The trian- 
gle in the lower panel is from llwata et al .1 j2009h. and the u pper 
limits at z = 1 .3 an d z ~ 4 come from ISiana et al. I j201Ch and 
IVanzella et al ] j201Ch respectively. 



and the intrinsic Lya emission, and there is currently no 
constraint on / osc at z ~ 0.3 — 2. Our results are consistent 
with the observations within the uncertainty. 

The evolution of the escape fraction of UV continuum 
at A rcst = 1300 — 1600 A, f^f, is shown in the middle panel 
of Figure [T] Not surprisingly, the general trend is similar to 
that of /esc™ j as /osc' is also determined by the dust con- 
tent and its distribution. However, since the path length 
of Lya photons before escaping can be longer than that of 
the UV continuum photons, so the absolute value of /^™ 
is expected to be smaller than that of ff s Y ■ For compari- 
son with observations, this figure also shows derived 
by using observed E (B — V) with Calzetti's extinction law 
jCalzetti et alj|200Ch . The attenuation by dust at 1600 A i s 
estimated by Ai 6 oo = 10 x E(B-V) (|Ouchi et al.ll2004h . 
Our result agrees well with the observations, which shows a 
moderate increase with redshift from z ~ 3 to z ~ 6. 

On the other hand, the escape fraction of the ionizing 



photons, f e °c, shows a different evolution, as shown in the 



lower panel of Figure [T] Unlike / 



esc 



and /esc i 



the does 



not change significantly with redshift, possibly due to ab- 
sorption of the ionizing photons by both neutral hydrogen 
and dust. In general, the mean gas density in galaxies in- 
creases with redshift, nu oc (1 + z) 3 , and the recombination 
rate is proportional to n|. Hence at high redshift, it be- 
comes difficult for the stellar radiation to ionize the gas. In 
addition, the mass ratio of gas to stars also increases with 
redshift, therefore, the ionizing photons can be more effec- 
tively absorbed by interstellar hydrogen gas at higher red- 
shift, leading to a low value of /'". On the other hand, at 
low redshift, although the mean gas density decreases, and 
galaxies are highly ionized by stars and UVB radiation, but 
the interstellar dust effectively absorbs the ionizing photons. 
As a result, /l™ does not increase at low redshift. 

The /esc 1 is not well constrained observationally, due 
to a number of uncertainties in the estimation, which in- 
clude the intrinsic SED, dust attenuation of the UV con- 
tinuum, and the IGM attenuation. In addition, most of ob- 
served objects were fainter than the detection limit. The /'™ 
from ou r simulations is clo se to the mean observed value at 
z ~ 3 (llwata et all [200$ ) and the upper limit at z ~ 4 
jVanzella et al]|201Ch. but somewh at higher than the upper 
limit at z = 1.3 (|Siana et aljboiCt ). This is likely due to the 
fact that our mod el galaxies are fainter than the sample of 
ISiana et all (|201Ch . 

Overall, our results of the escape fraction of the Lya 
and continuum photons are in general agreement with cur- 
rent observations. However, we note that the galaxy sample 
in this work is very small, and it is limited to progenitors 
of a MW-like galaxy. We will investigate the Lya and con- 
tinuum property of more galaxy populations from uniform 
simulations in larger volumes in future work. 



4.2 Relationship between the escape fraction of 
Lya and continuum photons 

In order to understand the different escape processes of 
Lya and continuum photons, we study in detail the rela- 
tion among the escape fractions of f^ a , /£c , and /i°" in 
Figure [2] We assess the strength of the correlation by fit- 
ting a simple power-law function, for instance, log/J^ = 
alog/egc" + /3, for galaxies at each redshift. The fitting is 
done on the mean value of each data point, with a bin size 
of 0.25 dex for the / e L sc Q . 

As shown in the figure, the scatter in the ff s Y — f^Z" 
relation (upper panel) appears to be smaller than that of the 
/esc 1 - /esc™ (lower panel). The slope of the former relation 
falls in the range of ~ 0.5 — 0.8 and does not evolve with 
redshift, while that of the latter decreases with redshift from 
a = 1.2 at z — to a = 0.1 at z = 7.2. This result suggests 
that /isc™ is more closely correlated with f^££ than with 
/esc 1 - The cause of the different relations is the difference in 
the escape processes. As mentioned in the previous section, 
both fe a Y and f^ a depend strongly on the dust content as 
it is the main absorber of Lya and UV photons, but ionizing 
photons can also be absorbed by neutral hydrogen gas, so 
/esc 1 depends on not only the dust, but the gas as well. 

To further demonstrate the escape processes of the Lya 
and UV continuum and their relation, it is useful to examine 
the "q" value, the relative strength of optical depth between 
Lya and UV continuum, q = t Lya It uv • Naively, one may 
think that since the traveling paths of Lya photons are much 
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Figure 2. The relation of the escape fractions between Lya and 
UV continuum at 1300 ^ A ^ 1600 A in the rest frame ( Up- 
per panel), and between Lya and ionizing photons (Lower panel) 
at different redshift. Solid lines arc absolute-lcast-dcviation fit- 
tings with a power-law function, log/^y (or fl°™ ) = alog/^c" + ft- 
The color of the points and lines corresponds to different redshift 
as specified in the color bar. 



longer than UV continuum photons due to numerous scat- 
tering processes, TL ya would be larger than tuv, resulting in 
high q values 1). Howeve r, recent observa t ions suggest 
that q ~ 1 - 2 at z = 0.3 - 3 jAtek et alj|2009l ; lHaves et afl 



l20ld ; lKornei et alllioiol ; iBlanc et alj|201lf ). Figure [3] shows 
the resulting "q" value from our calculations at different 
redshift. It has a mild variation over a wide redshift range 
2 = — 10, with a median value q ~ 1 — 2. This is in 



Figure 3. The evolution of ratio of optical depths between Lya 
and UV continuum (q = iXyc./'nyv') with redshift. The black open 
circles represent individual galaxies from the simulated sample, 
red filled circles are the corresponding median values at each red- 
shift, with the error bars showing the quartilcs. 



good agreement with observations, suggesting that photon 
escape is a complicated process affected by many factors 
other than scattering. At high redshift, the dust content 
is low, and it distributes compactly around galaxy center. 
Hence, although Lya photons can experience many scatter- 
ing processes over extended region, they can escape without 
dust absorption, leading to q ~ 1. At low redshift, since the 
amount of dust and metallicity increase, Lya photons are 
more effectively absorbed by dust than the UV continuum, 
resulting in somewhat higher q. 

4.3 Dependence of f csc on galaxy properties 

The relations between photon escaping process and galaxy 
properties and their evolution are not well known. Here we 
explore the dependence of photon escape fraction (/ c L s f , 
f^ s Y, and /es*) on various physical properties of the galaxies 
from our simulations. As shown in Figure 3J /<^c Q and /esc' 
appear to have similar relations with galaxy properties such 
as total mass, metallicity, dust mass, star formation rate 
(SFR), specific SFR (SSFR), and star formation efficiency. 
There is no strong correlation between photon escaping and 
the physical properties of the galaxy, although f^Z a an d /So 
show weak declining trend with the galaxy mass, metallicity 
and dust mass in the interstellar medium, albeit with large 
dispersions. Such a trend is in broad agreement with obser- 
vations that meta l enriched galaxies tend to have lower J^" 
(|Atek et al.ll2009l '). 

In galaxies with high SFR, there is a large supply of dust 
and neutral hydrogen gas for star formation, which can lead 
to efficient scattering and absorption of Lya photons. This 
trend is good agreement with a recent semi-analytic work of 
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Figure 4. Relation between photon escape fraction, f e / c a (upper panels), f^ s Y (middle panels), and fl^c (lower panels), and galaxy 
properties, including total galaxy mass (Mtot), gas mctallicity Z, dust mass (M(j ust ), star formation rate (SFR), specific SFR (SSFR), 
and star formation efficiency (csf) defined as the mass fraction of stars to gas. The color indicates different rcdshift as shown in the 
color bar. 



iGarel et all (|2012r i, which shows suppression of Lya escap- 
ing by dust in galaxies with high SFR. In addition, the Lya 
photons can frequently encounter dust in high-metallicity 
regions. Therefore, the Lya in galaxies of high SFR, metal- 
licity and dust mass can be effectively absorbed by dust, re- 
sulting in a low feZ a - Moreover, f^Z™ appears to roughly de- 
crease with halo mass. For example, at z = 3, it shows ~ 0.32 
at M h ~ 10 10 M Q to ~ 0.25 at M h ~ 10 " M w . This is 
good agreement with previous simulations bv lLaursen et all 
(|2009l ). in which they calculated fcZ a of nine galaxies at 
z — 3.6 and found similar relationship. 

However, /<!°", does not show any clear dependence on 
any of these properties. This again may be the result of a 
more complicated radiative process of the ionizing photons 
in dust and neutral hydrogen gas, and it may explain the 
lack of evolution with redshift seen in previous section. We 
note that some previous stu dies reported some correlation 
between /*°" and halo mass (R azoumov fc Sommer-Larsenl 
l20ld : lYaiima et ai1l2009l . l201lh . The discrepancy may be due 
to the different galaxy population in these studies. In our 
simulation, we focus on a MW-like galaxy, and our galaxy 
sample are mainly the pro genitors of such a galaxy. How- 
ever, in lYaiima et alj (|201ll ). their galaxy sample come from 
a uniform volume, which may represent the general popula- 
tion seen in the local universe. We plan to extend this study 
to a larger population of galaxies from larger, uniform sim- 
ulations. 
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Figure 5. Emissivity of ionizing photons as a function of emer- 
gent Lya luminosity. 



4.4 Contribution from LAEs to cosmic 
reionization 

The escaping of photons from early galaxies has a significant 
role in the reionization of neutral hydrogen. However, the 



estimation of ionizing photon emissivity (JVi on ) of observed 
LAEs has been difficult due to the ambiguity of the escape 
fraction of Lya and ionizing photons. Although the intrin- 
sic emissivity of ionizing photons is roughly related to SFR 
which can be derived from Lya luminosity, the estimation 
of the intrinsic Lya luminosity depends strongly on fcZ a - 
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Figure 6. Ionizing photon omissivity of LAEs at different red- 
shifts (colored symbols), in comparison with theoretical predic- 
tions of the number of ionizing photons needed for ionizing the 
IGM (solid lines). The emissivity is derived from the integration 
of the Lyo luminosity function converting from L^ ya to ionizing 
photon number. D ifferent color circle indicates different luminos- 
ity fu nction u sed: ICo wie et af| j2 010„ orange), Dcha.rveng et ah! 
(120081 . cy an). iGuaita et all <20ld . pur ple). iBIanc et al.l feoilT 
magenta), iGronwall et all J20Q7I. blue). lOuchi et all J2008l . red), 
lOuchi et all l|2010l . brown) and Iffibon et all l|20lT green). The 
filled and open circles represent the integration of the luminos- 



ity function to the detection limit and L^ y 



10 ergs s 



respectively. Triangles at 2 = 7 are the emissivity by using the 
luminosity function at 2 = 5.7 in lOuchi et all l|2008l ). The number 
of ionizing photons needed for ionizing the i nter-galactic medium 
(IGM) are calculated based on the model of lMadau et alj l|l999l ). 
N ion (s" 1 Mpc~ 3 ) = 10 47 4 C(1 + z) 3 , where C is the dumpiness 
of IGM. 



In addition, the emissivity of ionizing photons of galaxies 
can be greatly changed by / e °,?. Therefore, an accurate es- 
timate of feZ a and /eg" is critical to understand the cosmic 
reionization history. Here we directly estimate the Ai on of 
observed LAEs by using our results of the escape fractions 
of Lya and ionizing photons. 

We first evaluate the relation between AW and emer- 
gent Lhya- As shown in Figure [S] although Ai on roughly 
increases with Ll V ci , there is a large scatter. We fit the data 
from each redshift with a power-law function, logiVion = 
alogLLya + P- The slope a changes with redshift, ~ 
0.7 (z = 0.0), 0.8 (z = 3.1), 1.3 (z = 6.2) and 1.9 (z = 10.2). 
At higher redshift, typical galaxy mass decreases, and they 
have smaller /<!°", as seen in Figure[T] Moreover, the LL ya of 
these low-mass, high-redshift g alaxies can be largely boosted 
by the excitation Lya cooling (|Yaiima et al.|[2012bl ). There- 
fore, the slopes at higher redshift have higher values than 
their lower redshift counterparts. Next, we integrate the lu- 
minosity function (LF) of observed LAEs, and by converting 
Lhya to iVi on using the above fitting formula, we derive the 
emissivity of ionizing photons from these LAEs. 

Figure [6] shows the resulting emissivity at different red- 
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Figure 7. The stacked emissivities of ionizing photons normal- 
ized to total emissivity. Different color indicates different observa- 
tions used in the derivation of ionizing photo n emissi vity from the 



tions used in tnc derivation ot ionizing photon emissivity trom trie 
Lya lu minosity function: ICowie et all l|2010l . orange) , IBIanc et all 
bOlll . brown), [piichi et all fcuOSt red : 2 = 3.1, cy an : 2 = 3.7 
and p urple : 2 = 5.7). lpTidTlet**a*ll l|201Cll . blue) and lHibon et all 
i2012l , green). 



shifts, in comparison with theoretical predictions of the 
number of ionizing photons needed for ionizing the inter- 
galactic medium (IGM). Filled and open symbols indicate 
the emissivity derived from the integration of the LF up to 
the observational limit and L_L yQ = 10 39 ergs s _1 , which is 
fainter than detection limit of current observations by or- 
der 2-3 (e.g., ~ 10 42 ergs s" 1 : ICiardullo et~aHl2012D . Solid 
lines are the number of ionizing photons needed for ioniz- 
ing the inter- g alacti c medium (IGM) based on the model of 
iMadau et all (|l999l ). N ion (s^ 1 Mpc" 3 ) = 10 47 ' 4 C(1 + z) 3 , 
where C is the dumpi ness of IGM. Recent simulations sug - 
gest C ~ 3 - 10 (e.g., Illiev et all [20071 ; iPawlik et ai] 12009 ). 
Note that C — 1 corresponds to a homogeneous IGM, hence 
the dumpiness must be larger than unity in the cosmologi- 
cal model based on hierarchical formation. If the emissivity 
from LAEs is on or above the solid lines, then it means that 
LAEs alone can ionize the IGM. 

As shown in the figure, the emissivity is larger than the 
C = 10 line at z ~ — 4, and hence they can maintain 
the ionization of the IGM. At 2 = 5.7, the emissivity is 
marginally close to the C = 3 line. At 2 > 6, the emissivity is 
below the C — 1 line, and it is below the C — 3 line even if we 
integrate the LF down to LL ya — 0. However, the observed 
LF at z > 6 could be suppressed by IGM scattering. Hence, 
the LF at z ~ 5.7 is frequently used as an intrinsic LF at 2 > 
6 (e.g.. iKashikawa et alllioill ). motivate d from no observe d 
evolution of the LF from z ~ 3.1 to 5.7 (|Ouchi et al.ll200&t ). 
So here we esti mate the emissivity at z = 7 by using the 
LF at 2 = 5.7 of lOuchi et all |2008l ). However, the derived 
emissivity is below the C — 3 line even considering the upper 
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limit. Therefore, we suggest that LAEs alone cannot ionize 
the IGM at z > 6. 

It has been suggested that other populations (e.g., 
LBGs, Pop III stars, faint quasars) may be needed to ionize 
the neutral hydrogen at z > 6. It was shown t hat the emis- 
sivity from LB Gs can be larger t han LAEs (|Yaiima et alj 
12003 ). However. lOuchi et alj <|2009f > showed that the emissiv- 
ity from LBGs at z ~ 7 falls below the C = 1 line. The other 
possibility to ionize the IGM by LAEs or LBGs is a correc- 
tion of the fai nt-end slope of LF. In observations of L AEs 
at z > 3 (e.g.. lOuchi et al.ll2010l ; ICiardullo et alj|2012» . the 
slope is assumed to be —1.36 ~ —1.65, w hich is derived 
from recent deep survey at lower redshift (e.g. JCassata et all 
l201lT ). because observations of high redshifts can only trace 
the bright-end of the LF. If this slope is steeper (e.g., 
a < —2), then the contribution from faint galaxies be - 
comes large (e.g.. llBouwens et al.ll20l"ll ; Ijaacks et al.ll201ll ). 
giving rise to higher emissivity which can ionize the IGM. 
The next generation telescopes which have higher sensitivity 
(e.g., JWST, GMT, TMT) will be able to measure the slope 
more precisely. 

Figure [7] shows the stacked emissivity as a function 
of LLya from observations of LAEs. The emissivity in- 
creases with redshift, as the slope of fitting function be- 
tween iVi on and LLya becomes steeper at higher redshift. 
About half of the ionizing photons comes from LAEs of 
4.2 x 10 41 sg L Lv g s£ 4.2 x 10 42 ergs s 
lOuchi et alj|2010h. 2.4 x IP 4 * $ L j.„„ < 2.4 x 
at z = 7 (LF by Hibon et al.ll2012j ). Hence, typical observed 
LAEs around L* Lya can be the main contributors to the 
IGM ionization among the LAE population. On the other 
hand, at z — 3, half of ionizing photons come from LAEs of 



at z = 6.6 (LF by 
10 42 ergs s" 1 



1.3 x 10 40 < Llvc < 7.5 x 10 41 ergs s" 1 (LF by lOuchi et all 



l200ct ). Therefore, LAEs fainter than dete ction limit of most 
observation (e.g., LL ya ~ 10 42 ergs s _1 : lOuchi et aI.ll2008T ) 
may contribute significantly to the IGM ionization. 



5 DISCUSSION 

5.1 Star formation rate from emergent Lya and 
UV flux 

The observed Lya and UV continuum fluxes are frequently 
used to derive SFR. However, such method is subject to 
the uncertainty of the photon escape fraction, and depend- 
ing on the model, the derived SFR may differ by orders of 
magnitude. To illustrate this point, we show in Figure[8]the 
comparison of SFRs derived from Lya and UV continuum 
flux at different redshift. 

The SFRs of some bright LAEs at z = 3.1 from our 
simulations are consi stent with the recent observation of 
ICiardullo et al.1 (|2012i ). but the dispersion is large. Since 
is higher than f^Z", the derived SFR UV is somewhat larger 
than SFR Lya , and close to the true SFR. Most of our model 
galaxies are fainter than the observational sample, so our re- 
sults represent the SFR relation of faint LAEs, and suggest 
that it roughly can be fitted by the same power-law relation. 

On the other hand, at z > 6, SFR Lyct becomes larger 
than SFR UV . At this high redshift, the contribution of exci- 
tation Ly a cooling to the emerg ent Lya luminosity becomes 
dominant (|Yaiima et al ] |2012bf ). The Lya cooling due to ex- 
citation is no longer directly related to SFR, so it can result 




-3 -2-10 1 
log SFR(Lya) 



Figure 8. Comparison of SFR derived from emergent UV con- 
tinuum and Lya photons. The esti mation of SFR fro m UV con- 
tinuum is based on Equation (1) in iKennicut j lll998h. The SFR 
from Lya is estimated using equation (2) in iKennicuttJ il998T ) 
with the assumption of case B. The cross symbols are the de- 
rived SFRs of spectro scopically confirmed LAEs at z = 3.1 from 
ICiardullo et al The color indicates different redshift as 

shown in the color bar. 



in overestimate of SFR, because the derivation of SFR by 
Lya flux as in iKennicuttl (|l99cf ) considers only stellar ra- 
diation. Deep multi-band surveys by the next generation 
telescopes (e.g., GMT, TMT) may trace such a tail of high 
SFR by Lya which may be evidence of dominant Lya source 
by excitation cooling process. Moreover, the contribution 
from Lya cooling due to excitation becomes larger in lower- 
mass galaxies due to lower star formation rate. Since SFR 
roughly increases with galaxy mass, as a result, the gradient 
between SFR UV and SFR LyQ becomes large. If we fit them 
by a power-law function, the power-law index is a — 0.6 
(z = 3.1), and 1.2 (2 = 10.4) 

5.2 Dependence of ISM dumpiness 

It was suggested by iNeufeldl (|l99lf ) that large EW of Lya 
line may be produced by a clumpy interstellar media, be- 
cause Lya photons can be scattered by neutral hydrogen 
on the surface of clouds before dust absorption, and they 
can escape without much dust absorption in the clumpy 
dusty media. On the other hand, UV continuum photons 
can be absorbed by dust passing through the clouds. As 
a result, f^Z a can b e larger than f^, causing higher EW. 
lHansen fc Ohl (|2006T ) showed that Lya can effectively escape 
in the idealized clumpy medium by RT simulations. 

Our simulations provide a unique opportunity to study 
the effect of ISM dumpiness on the escape of Lya and UV 
continuum photons. Figure [5] shows the ratio between f^Z™ 
and /e^J as a function of the dumpiness of the interstellar 
gas, defined as C =< p 2 > / < p > 2 . We find that the / e L a r 
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by combining our simulations with the observed luminosity 
functions of LAEs at different redshift. We found that ioniz- 
ing photons from LAEs can maintain the ionization of IGM 
at z ~ — 4, but they are not sufficient to ionize IGM at 
2 > 6, other populations (e.g., LBGs, Pop III stars, faint 
quasars) may be needed to ionize the neutral hydrogen. 

Moreover, we found that the star formation rates de- 
rived from emergent Lya and UV may subject to uncer- 
tainty of orders of magnitude depending on the photon es- 
cape fractions. This would have significant implications in 
the interpretation of observational data. Furthermore, we 
found no correlation between the EW of Lya line and the 
clumpiness of the ISM suggested by previous studies. 

Overall, our results of the escape fraction of the Lya and 
continuum photons are in general agreement with current 
observations. However, we note that the galaxy sample in 
this work is very small, and it is limited to progenitors of a 
MW-like galaxy. We will investigate the properties of Lya 
and continuum photons of a larger galaxy populations from 
uniform simulations in larger volumes in future work. 



Figure 9. The ratio of escape fraction of Lya to UV continuum 
photons as a function of the clumpiness of the interstellar gas, 
which is defined by < p 2 > / < p > 2 , where p is gas density. The 
color indicates different redshift as shown in the color bar. 



not correlate with the clumpiness. In our simulations, the 
Lya photons can be created in high-density clouds, because 
star formation can occur in such a region and ionize the gas 
partially. In such a situation, the Lya photons are effectively 
absorbed by dust at the origin place. Hence, the cannot 
be higher than /JJj even in clumpy structure. Our results 
suggest that there is no correlation between the Lya EW 
and the clumpiness of the ISM. 



6 SUMMARY 

To summarize, we have investigated the escape of Lya, 
non-ionizing UV continuum, and ionizing photons from 
galaxies by combining cosmological SPH simulations and 
three-dimensional multi-wavelength radiative transfer calcu- 
lations. We found that the escape process differs significantly 
from one photon species to another. correlates more 



strongly with /^J than with /i™. Both /, 



Lya 



and fcic have 



similar evolution trend with redshift, and they show simi- 
lar, weak dependence on galaxy properties such as galaxy 
mass, metallicity, dust content, star formation rate, specific 
star formation rate, and star formation efficiency. However, 
fl°c does not show evolution with redshift, and it shows 
no clear correlation with galaxy properties. These different 
behaviors may be explained by the different escaping mech- 
anisms. Both and depend strongly on the dust 
content as it is the main absorber of Lya and UV photons, 
but ionizing photons can also be absorbed by neutral hydro- 
gen gas, so /cs" depends on not only the dust, but the gas 
as well. 

In addition, we estimated the ionizing photon emissiv- 
ity of LAEs and their contribution to the ionization of IGM, 
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